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We show that all of the relevant features of a phase transition can be determined using a non 
order parameter field which is a physical state of the theory. This fact allows us to understand the 
deconfinmg transition of the pure Yang-Mills theory via the physical excitations rather than using 
the Polyakov loop. 
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INTRODUCTION 

Some of the relevant intrinsic properties of the phase 
transition are best investigated using the order param- 
eter of the theory. Sometimes the order parameter can 
be a non physical object or not experimentally accessi- 
ble. The pure Yang-Mills deconfinement phase transition 
constitutes a time honored example. Here the order pa- 
rameter is the Polyakov loop which is not a physical state 
of the theory while the physical field is the glueball field. 
We argue that it is always possible to extract all of the 
relevant information about a specific phase transition us- 
ing a non order parameter field which is a physical state 
of the underlying theory. 

We predict the universal behavior of this field near and 
at the phase transition. More specifically we predict a fi- 
nite drop in the mass of the singlet field at the phase 
transition. Our predictions seem to reproduce surpris- 
ingly well the general features of some lattice data and 
can be put to further and more stringent tests via new 
first principle lattice simulations. Finally, we define new 
physical quantities able to characterize the phase transi- 
tion. 



A RENORMALIZABLE THEORY AND 
CRITICAL BEHAVIOR 

We consider a temperature regime close to the phase 
transition. In order for our results to be as universal as 
possible we study a renormalizable Lagrangian contain- 
ing a field neutral under the global symmetries and the 
order parameter field of theory as well as their interac- 
tions. The protagonists of our theory are two canonically 
normalized fields h and x- The field h is a scalar of Zn 
while x transforms according to x ~ > z X with z € Zn- 

Here we consider the case N = 2. This theory is suit- 
able for understanding the deconfining phase transition 
of 2 color Yang-Mills which has been heavily studied via 
lattice simulations 1,2]. A general, renormalizable, po- 
tential is: 
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The coefficients are all real with go — 0. Stability requires 
A > and g± > 0. It is straightforward to generalize this 
potential to Zn- Since we will relate x to the Polyakov 
loop, we take it to depend only on space coordinates. 
This leads to a mixed situation in which the physical 
states depend also on time. We will also briefly comment 
on what happens when the order parameter is taken to 
be four dimensional. The Lagrangian we use to define 
our Feynman rules is: 
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We start our investigation at temperatures below T c . 
Since the phase transition is of second order, m x vanishes 
at the transition point. All of the coupling constants are 
such that h and x have zero vacuum expectation value 
in this regime. 

We focus our attention on the h self energy at the one 
loop level at a given non zero temperature: 
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The double lines indicate x fields while single lines stand 
for h ones. We are not considering tadpole diagrams 
which we require to vanish at each given loop order by 
adding a counterterm linear in h. The number in front 
of each diagram is the associated combinatorial factor. 

We consider the limit in which the h mass is much 
larger than the temperatures in play, hence all of the 
loops containing only h fields are then infrared finite at 
any temperature. The second diagram in eq. is, how- 
ever, linearly ultraviolet divergent but this divergence is 
absorbed in the respective mass counterterm. After all 
the ultraviolet divergences have been taken into account, 
we can show that the finite temperature corrections due 
to h are Boltzman suppressed. We are left with the fol- 
lowing infrared divergent graph for T — > T c , 



-o- 



T ■ 



9\ 



167rmT, 



(4) 



and this diagram constitutes the relevant one loop cor- 
rection to the h mass which reads: 
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This is the way the nearby phase transition is directly 
felt by the non order parameter field, ft shows that the 
mass of the singlet field must decrease fast close to the 
phase transition. An unpleasant consequence is that this 
one loop result breaks down at the transition point due 
to the infrared singularity. 

Besides, since h is not the order parameter, its cor- 
relation length (i.e. 1/m) is not expected to diverge at 
the phase transition. To cure the infrared behavior we 
need to go beyond the one-loop approximation and con- 
sistently resum the following set of diagrams: 
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Each diagram is infrared divergent but when re- 
summed to is finite at T and reads: 
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At T = T r we have 



2/7rt\ 2 9l 
TO [I c) = TO — 



(6) 



(7) 



Hence we predict that close to the phase transition the 
singlet state has a decreasing mass. The drop at the 
phase transition point is the ratio between the square 
of the coupling constant governing the interaction of the 
singlet state with the order parameter (gi ) and the order 
parameter self interaction coupling constant (A). 
The slope is: 
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with Am 2 = to 2 (T c ) — to 2 = g 2 /X. This slope encodes 
the critical behavior of the theory even though it is con- 
structed using the non order parameter field. For ex- 
ample if to 2 vanishes as (T c — T) v close to the phase 
transition, then T>~ scales with the exponent {v/2 — 1). 
When considering an 0{N) symmetry group rather than 
Z<i the resummation procedure becomes exact in the large 
N limit @. 

When T > T c x develops a vacuum expectation value 
which will induce one also for h: 
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where we denoted the mass of \ m the broken phase 
by M 2 = 2|to x | 2 and considered small h fluctuations 
(i.e. we kept only the gi interaction term for h). This 
is in agreement with the results found in The fields 
X and h mix in this phase and the mixing angle 9 is 
proportional to g 2 /m 2 . The mixing can be neglected 
within the present approximations. Like for T < T c , we 
consider also now only the effects due to x loops for the 



h propagator. Due to symmetry breaking we now have 
the trilinear x couplings: 



A 
3! 



vx 



(10) 



which substantially affect the analysis and the results in 
this phase. At the one loop level the diagram to compute 
is again the one in eq. (@J with m x replaced by M x . So 
we predict a drop on the right hand side of T c while 
this diagram is clearly infrared divergent. Curing the 
divergence now is more involved due to the appearance 
of the trilinear x coupling. A typical diagram in the 
consistent set of diagrams we consider is of the form: 
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We choose this subset of diagrams since the associate in- 
finite sum can be performed exactly 0| and it has knowl- 
edge of the onset of chiral symmetry breaking via the 
trilinear vertices. 
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The infrared divergence has been cured for T > T c . We 
also see that m(T c ) from the right hand side of the transi- 
tion equals exactly the one from the unbroken side of the 
transition. This remarkable result does not hold order 
by order in the loop expansion but only when an infinite 
sum of the diagrams is performed. The h mass square is a 
continuous function throughout the phase transition and 
the associated correlation length remains finite. However 
the slope of the h mass near the phase transition from 
the right hand side does not mirror the one on the left. 
This fact is due to the onset of spontaneous symmetry 
breaking communicated to h via the trilinear x interac- 
tion term. Indeed: 
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and we deduce the relation: 
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The h mass drops on the left hand side of the phase 
transition and rises on the right one. The slope on the 
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FIG. 1: Left panel: behavior of the mass of the singlet field 
close to the phase transition as function of the temperature. 
Right panel: Lattice data; figure from |(J. 



right side (broken phase) is less singular with respect to 
the one on the left. More specifically the scaling expo- 
nent for T> + is (y — 1). Deviations from such a behavior 
directly measure departures from the second order char- 
acter of the phase transition. In the figure we schemati- 
cally represent the behavior of the h mass as a function 
of the temperature in units of the critical temperature 
for ex (Tr — T). The similarity of our result with the 
lattice data Q (right panel of Fig^l is surprising. Re- 
call, though, that our results are directly applicable to 
the two color Yang-Mills theory, while the data in 6] are 
obtained for the three color Yang-Mills theory This re- 
semblance is in agreement with the almost second order 
character of the three color deconfining phase transition. 
Another difference is that the mass in seems to be 
the screening one while we are plotting here the physical 
mass. 

It is of common practice to isolate the order parame- 
ter and more generally the light degrees of freedom since 
they are expected to be the relevant degrees at low en- 
ergies. While this is certainly correct we have explicitly 
shown how all of the essential features of the phase tran- 
sition are encoded in the non order parameter fields of 
the theory as well. Our analysis is useful whenever the 
order parameter is neither a physical quantity nor phe- 
nomcnologically accessible (e.g. the Polyakov Loop). 

We also studied what happens if the order parame- 
ter depends on the time coordinate. The h mass has the 
same behavior near the transition as in the previous case. 
This can be easily understood since the infrared physics 
is dominated by the zero mode of %■ The major differ- 
ence with the previous analysis is that a physical width 
develops describing the decay of h into two % fields at 
non zero temperature 0. 



DECONFINEMENT AND CONCLUSIONS 

The SU(N) deconfmement phase transition is a hard 
problem. Importance sampling lattice simulations pro- 
vide crucial information about the nature of the temper- 
ature driven phase transition for 2 and 3 colors Yang- 
Mills theories with and without matter fields (see [7J for 



the features of this phase transition. At zero tempera- 
ture SU(N) Yang-Mills theory is asymptotically free and 
the physical spectrum of the theory consists of a tower 
of hadronic states referred to as glueballs and pseudo- 
scalar glueballs. The theory develops a mass gap and the 
lightest glueball has a mass of the order of few times the 
confining scale. 

At nonzero temperature the Zn center of SU (N) is a 
relevant global symmetry |20|. and it is possible to con- 
struct a number of gauge invariant operators charged 
under Zn among which the most notable one is the 
Polyakov loop: 
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V denotes path ordering, g is the SU(N) coupling con- 
stant and x is the coordinate for the three spatial di- 
mensions while t is the euclidean time. The £ field is 
real for N — 2 while otherwise complex. This object 
is charged with respect to the center Z N of the SU(N) 
gauge group |20| | under which it transforms as t — > z£ 
with z € Zn- A relevant feature of the Polyakov loop 
is that its expectation value vanishes in the low temper- 
ature regime while is non zero in the high temperature 
phase. The Polyakov loop is a suitable order parameter 
for the Yang-Mills temperature driven phase transition 

E3 

Pisarski used this feature [13j to model the Yang-Mills 
phase transition via a mean field theory of Polyakov 
loops. The model is referred to as the Polyakov Loop 
Model (PLM). Recently some interesting phenomenolog- 
ical PLM inspired models aimed to understand RHIC 
physics were constructed 0, |52 • 

Here we consider pure gluon dynamics. This allows us 
to have a well defined framework where the Zn symmetry 
is exact. The hadronic states of the Yang-Mills theory are 
the glueballs. At zero temperature the Yang-Mills trace 
anomaly has been used to constrain the potential of the 
lightest glueball state H 0: 



V[H] = I L In 



H 
A 7 



(17) 



A is chosen to be the confining scale of the theory and 
H is a mass dimension four field. This potential cor- 
rectly saturates the trace anomaly when H is assumed 
to be proportional to Tr [G^G^] and G^ v is the stan- 
dard Yang-Mills field strength. 

In 0] a concrete model was proposed able to trans- 
fer the information about the Yang-Mills phase transi- 
tion encoded in the Zn global symmetry to the hadronic 
states of the theory. This model is constructed using 
trace anomaly and the Zn symmetry. The model in is 
supported by recent theoretical investigations ji^] . Once 



4 



the relation between the fields H and h, and £ and \ 
is made the relations between the couplings of the La- 
grangian for H and £ constructed in [4J and the renor- 
malizable one presented here can be obtained. More 
specifically, our renormalizable theory is a truncated (up 
to fourth order in the fields) version of the full glueball 
theory. For example one can take the following relation 
between H and the glueball field h: 

E -w( 1+ 7m*)- (18 » 

Here (H) = A 4 /e is the vacuum expectation value of 
the glueball field below the critical temperature and c is 
a positive dimensionless constant fixed by the mass of 
the glueball. For the £ field we have x = V^^ with k 
a mass dimension two constant (at high temperature is 
proportional to T 2 ). 

Our present results are simply the higher loop correc- 
tions to the glueball model presented in jj] and can be 
immediately applied to the two color Yang-Mills phase 
transition. We confronted already our theoretical results 
with lattice computations of the glueball mass behavior 
close to the phase transition studied in Q for three colors. 
Our analysis not only is in agreement with the numerical 
analysis but allows us to provide a better understanding 
of the physics in play. More recent lattice results should 
provide a still more detailed test of our predictions. 

Finally, using a renormalizable theory containing a sin- 
glet heavy field uncharged under a global symmetry and 
the associated order parameter with their possible inter- 
actions we demonstrated that the information about the 
phase transition can be extracted via the singlet field. 
The transfer of information from the order parameter to 
the heavy field is in fact complete. Moreover, if we con- 
sider the large N limit for the O(N) theory the analysis in 
the unbroken phase becomes exact. Thanks to the result 
established here and envisioned in Q we can use directly 
the knowledge about h, or more specifically of some of its 
properties defined in the text, to characterize the phase 
transition. We considered a Z2 symmetry group appro- 
priate for example to discuss the deconfining phase tran- 
sition for 2 color Yang-Mills theory. The generalization 
to a general group of symmetries is straightforward. 
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